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Abstract 

We measured the trigonometric annual parallax of H2O mascr source associated with the massive star- 
forming regions IRAS 06061+2151 with VERA. The annual parallax of 0.496 ±0.031 mas corresponding 
to a distance of 2.02^^2 kpc was obtained by 10 epochs of observations from 2007 October to 2009 
September. This distance was obtained with a higher accuracy than the photometric distance previously 
measured, and places IRAS 06061+2151 in the Perseus spiral arm. We found that IRAS 06061+2151 also 
has a peculiar motion of larger than 15 km s _1 counter to the Galactic rotation. That is similar to five 
sources in the Perseus spiral arm, whose parallaxes and proper motions have already been measured with 
higher accuracy. Moreover, these sources move at on average 27 km s _1 toward the Galactic center and 
counter to the Galactic rotation. 

Key words: astrometry — masers (H2O) — ISM: Hll regions — techniques: high angular resolution 
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1. Introduction 

The annual parallax is the most direct way to accu- 
rately measure the distance to a target source without any 
assumption. Very Long Baseline Interferometer (VLBI) 
provides the highest accuracy of position determination 
among astronomical instruments. It is possible to deter- 
mine the position of target sources with an accuracy of 
less than 1 mas with the VLBI observations. The phase- 
referencing VLBI has been utilized in recent astrometry to 
measure the annual parallaxes of target sources. For ex- 
ample, the results of trigonometric parallax measurements 
of H2O and CH3OH maser sources using the NRAO Very 
Long Baseline Array (VLBA) were reported by Xu ct al. 
(2006); Hachisuka et al. (2006, 2009); Reid et al. (2009a). 

The Japanese VLBI Exploration of Radio Astrometry 
(VERA) is an instrument dedicated to carrying out VLBI 
astrometry. It aims to figure out the Galactic 3D structure 
by measuring the trigonometric parallaxes and the proper 
motions of maser sources with the phase-referencing 
VLBI. To carry out effective phase-referencing observa- 
tions, the dual-beam system was developed (Kobayashi 
ct al. 2008). The results of the trigonometric paral- 
laxes of several H2O maser source (e.g. Orion-KL, IRAS 
06058+2138) with VERA have already been reported by 
Choi et al. (2008); Hirota et al. (2007, 2008a, 2008b); 



Honma et al. (2007); Imai et al. (2007); Nakagawa et al. 
(2008); Oh et al. (2010); Sato ct al. (2008, 2010). 

It is believed that our Galaxy is composed of several spi- 
ral arms. In order to understand the spiral structure and 
the rotational kinematics of the Galaxy, it is important to 
study the structure of each spiral arm. So, measuring the 
accurate distance to the astronomical objects and investi- 
gating their kinematics in the spiral arms are essential for 
studying spiral arm structure. The Perseus spiral arm is 
the nearest one to us outward from the Sun in the Galaxy 
and consists of many star-forming regions. Also, the dis- 
tances to six sources have previously been measured with 
high accuracy in previous studies (Moellenbrock et al. 
2009; Moscadelli et al. 2009; Oh et al. 2010; Reid ct al. 
2009a; Sato et al. 2008; Xu ct al. 2006). It turns out that 
some of these sources (e.g. W30H, IRAS 00420+5530) 
showed their peculiar motion with respect to the Galactic 
rotation. However, it is necessary to increase the number 
of samples because the current number is insufficient to 
indicate that their peculiar motions relate to the motion 
of the Perseus spiral arm. 

IRAS 06061+2151 is a massive star-forming region 
located in the Gemini OBI cloud complex. Several 
sources are associated with IRAS 06061+2151, for ex- 
ample, a radio continuum source associated with the 
ultra compact Hll region, G188. 794+1. 031 (Kurtz et 
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Table 1. VERA observations of IRAS 06061+2151 



Epoch ID 


Date 


Day of year Synthesized beam 




T 

sys 








Iriki 


Mizusawa Ogasawara Ishigakijima 






[mas x mas (°)] 




[K] 



1 


2007/Oct 23 


2007/296 


1.80 x 1.12 (128.5) 


150 


- 240 


110 - 


150 


320 


- 600 


240 - 


550 


2 


2007/Dec 15 


2007/349 


1.83 x 1.13 (127.4) 


110 


- 170 


110 - 


530 


150 


- 210 


250 - 


550 


3 


2008/Feb 15 


2008/046 


1.81 x 1.13 (128.9) 


100 


- 160 


100 - 


190 


150 - 


2600 


150 - 


300 


4 


2008/ Apr 20 


2008/111 


1.84 x 1.13 (127.4) 


130 


- 200 


130 - 


200 


160 


- 270 


250 - 


550 


5 


2008/ Jun 15 


2008/167 


1.79 x 1.12 (128.5) 


500 - 


12000 


180 - 


260 


300 


- 550 


120 - 


200 


6* 


2008/Sep 29 


2008/273 


1.74 x 1.09 (129.7) 






120 - 


200 


450 


- 850 


350 - 


1500 


7 


2008/Nov 02 


2008/307 


1.86 x 1.13 (125.7) 


250 - 


2300 


140 - 


280 


200 


- 380 


300 - 


750 


8 


2008/Dec 14 


2008/349 


1.90 x 1.13 (133.4) 


100 


- 140 


100 - 


130 


200 


- 380 


180 - 


1400 


9 


2009/Feb 15 


2009/046 


1.86 x 1.14 (130.7) 


150 


- 700 


110 - 


180 


240 - 


5000 


210- 


540 


10 


2009/May 11 


2009/131 


1.82 x 1.13 (128.5) 


100 


- 145 


200 - 


340 


150 


- 240 


200 - 


450 


11 


2009/Sep 07 


2009/250 


1.89 x 1.15 (128.5) 


190 


- 360 


250 - 


580 


290 


- 600 


210 - 


440 



*: Wc do not include the data of the epoch 2008/273 for the astromctric analysis (Sect. 2) 
**:In the column of T SJ 
in each epoch. 



; the observation at each station and 



al. 1994), a sub- millimeter wavelength source JCMTSF 
J060907.1+215037 (Di Francesco et al. 2008), a near 
infrared (NIR) cluster GL 5182, the C 18 line emission 
(Saito et al. 2007). Moreover, this region shows the 
H2O maser emission (Kompe et al. 1989). Thus the 
structure and the kinematics of a molecular gas around 
this region have been well studied with high spatial res- 
olution observations. Especially, Motogi et al. (2008) 
found the bipolar morphology and the expanding motion 
suggested by the relative proper motion measurement of 
H2O masers with the Japanese VLBI Network (JVN: the 
maximum angular resolution is about 1 mas at 22 GHz). 
The photometric distance to the Gemini OBI cloud com- 
plex was measured by the previous studies, and the dis- 
tance range was from 1.5 kpc to 3.5 kpc (Carpenter et 
al. 1995; Humphreys 1978; Moffat et al. 1979). IRAS 
06058+2138 whose parallax distance has already been 
measured as approximately 2 kpc by VLBI observation, is 
known to be located in the Perseus spiral arm and located 
near IRAS 06061+2151 on the celestial sphere. Therefore, 
although the uncertainty of the photometric distance is 
much larger, there is a possibility that IRAS 06061+2151 
may also be located in the Perseus spiral arm. 

In this paper, we report on the result of the annual 
parallax measurement of the H2O maser source in IRAS 
06061+2151 with VERA, and discuss the motion in the 
Galaxy. 

2. VLBI observations with VERA 

VERA observations of the H2O maser source (616 — 523, 
22235.080 MHz) associated with IRAS 06061+2151 have 
been carried out since 2007 October to 2009 September at 
2- to 3-month intervals. We present the results from a to- 
tal of 11 observing epochs: 2007/296, 2007/349, 2008/046, 
2008/111, 2008/167, 2008/273, 2008/307, 2008/349, 
2009/046, 2009/131, and 2009/250; hereafter an observ- 
ing epoch is denoted by year/day of the year (DOY). All 4 
VERA stations, providing a maximum baseline length of 



2270 km, participated in 10 epochs, but only 3 VERA sta- 
tions (Mizusawa, Ogasawara, and Ishigakijima: see Figure 
1 of Pctrov et al. (2007)) participated in the 2008/273 
epoch. Therefore, we do not include the data of the epoch 
2008/273 in order to conduct a more accurate astrometric 
analysis. 

Observations were made in the dual-beam mode and 
we observed simultaneously (Kawaguchi et al. 2000; 
Honma et al. 2008a) the H2O maser source in 
IRAS 06061+2151 (J2000.0 a = 06 h 09 m 06?97458±0?006, 
5 = +21°50'4i:'4045 ± 0."05; Migenes et al. 1999), 
and the reference source J0603+2159 (J2000.0 a = 
06 h 03 m 51?55709±0?00008, 5 = +21°59'37."6976±0f'0035; 
Fey et al. 2004) . The separation angle between those two 
was l.°22. J0603+2159 was detected with a peak flux den- 
sity ranging from about 60 to 140 mJy beam -1 throughout 
all epochs. Left-handed circular polarization was received 
and sampled with 2-bit quantization and filtered using 
the VERA digital filter unit (Iguchi et al. 2005). The 
data were recorded onto magnetic tapes at a rate of 1024 
Mbps, providing a total bandwidth of 256 MHz which was 
divided into 16 IF channels (16 MHz each). 8 MHz of the 
one IF was assigned to IRAS 06061+2151 and the other 
15 IFs were assigned to J0603+2159. A bright contin- 
uum source, DA193, was observed every 90 minutes for 
bandpass and delay calibration. The system noise tem- 
peratures (T sys ) including atmospheric attenuation were 
depending on the weather conditions and the elevation 
angle of the observed sources. We list the value of T^ 1 
as well as the typical synthesized beam sizes (major and 
minor axes) and their position angles at all epochs in Table 
1. The aperture efficiencies of the antennas ranged from 
45 to 52 %, depending on the stations (see, the VERA sta- 
tus report 2 ). A variation of the aperture efficiency of each 
antenna as a function of the elevation angle was confirmed 
to be less than 10 %, even at the lowest elevation in the 

1 We list the minimum and the maximum of T sys during the ob- 
servation at the each station and the each epoch in Table 1. 

2 http: / /veraserver. mtk.nao.ac.jp/restrictcd/CFP2009/status09.pdf 
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2008/307 (1.61) 




2009/131 (3.43) 



4 2 0-2 


-4 


4 2 


-2 -4 


4 2 -2 -4 


X [mas] 




X [mas] 




X [mas] 


10 20 30 40 


50 60 


10 20 30 40 


60 eo 70 


2 4- 6 8 10 12 




Peak intensity 




Peak intensity 


Peak intensity 




61.078 




74. 380 


12. 117 




[Jy/beam] 




[Jy/beam] 


[Jy/beam] 



Fig. 1. The phase-referenced map of the maser spot at i>LSR = —7.795 km s -1 in three epochs. Epochs when the maser spot 
was observed and the ratio of the maser spot size whose flux has more than 30% of brightest component in the map to beam size 
are shown on the top of each map, as 'DOY (spot size)'. The most right figure was not included in our result because spot size 
of the maser spot exceeds the threshold for conducting an astrometric analysis in the presented study (see Sect. 4.1, and 5.1). 



observations (~20°). Correlation processing was carried 
out on the Mitaka FX correlator (Chikada et al. 1991) 
on the NAOJ Mitaka campus. For H2O maser lines, the 
spectral resolution was set at 15.625 kHz, corresponding 
to a velocity resolution of 0.21 km s _1 . So, the bandwidth 
of 8 MHz corresponds to the velocity coverage of 107.5 km 

3. Data reduction with VEDA 

We conducted the data reduction with the software 
package VERA Data Analyzer (VEDA). VEDA has been 
developed for the astrometric analysis of data observed 
with the VERA dual-beam system by the software devel- 
opment group in Mizusawa VLBI observatory of NAOJ. 

At first we conducted the averaging of time and fre- 
quency for the data of the position reference source 
J0603+2159 in order to reduce the data size. After data 
averaging, we carried out the amplitude calibration and 
the bandpass calibration for the maser source in IRAS 
06061+2151 and J0603+2159 independently. The ampli- 
tude calibrations were performed based on the (T sys ) data 
monitored during the scan. The bandpass calibrations 
were performed utilizing auto-correlation spectra of the 
blank sky data and the bright calibrator source (DA193). 

The fringe search was conducted with the visibility data 
of the fringe finder DA193 in order to determine the 
delay and the delay-rate offset. Next the fringes were 
searched for J0603+2159 by utilizing the station solu- 
tion determined by DA193, and the image of J0603+2159 
was obtained by the deconvolution (CLEAN) and the self- 
calibration algorithm. 

The absolute position of the maser source was obtained 
by applying the phase solutions of J0603+2159 to the 
maser source (phase-referencing method). For the maser 
source, the Doppler tracking were made. Finally, we made 
images for each velocity channels, which were selected so 
that LSR velocities of maser features reported in Motogi 
et al. (2008) are fully covered, in order to search for maser 



spots. Each of the velocity channels was imaged with a 
field of view of 2" x 2" with 1024 grids x 1024 grids by 
CLEAN with a threshold for signal to noise ratio of 7 (7a 
detection limit). For each detected maser spot, the images 
of 10 mas x 10 mas were obtained with 512 grids x 512 
grids. 

After the imaging process, we conducted the astromet- 
ric analysis using the result of phase-referenced analysis. 
For the astrometric analysis, the maser spots without ex- 
tended structures should be selected so as to accurately 
determine the position. So, we can set the threshold for 
the structure of maser spot using the parameter "spot 
size" with VEDA. The spot size is defined as the ratio of 
an area whose flux is larger than 30 % of the peak flux 
density to the synthesized beam size. We set the limit 
of 5.0 for the spot size in imaging process, so the maser 
spots whose size extend beyond 5.0 times the beam size 
were rejected. As an example, we show channel maps of 
the same spot detected with a different spot size at several 
epochs in Figure 1. We discussed the details of the spot 
size in Sect. 5.1. 

4. Results 

The cross-power spectrum of the H2O masers to- 
ward IRAS 06061+2151 from 2007/296 to 2009/250 with 
Mizusawa-Ishigakijima baseline (2270 km) are shown in 
Figure 2. The component around Ulsr = — 8 km s _1 
showed the stable emission through the entire epoch. The 
systemic velocity {v sys ) of IRAS 06061+2151 is -1.6 + 0.2 
km s _1 , which is obtained by 12 CO observations of IRAS 
sources (Wouterloot et al. 1989). In total we detected 142 
maser spots whose LSR velocities are ranged from —19.2 
km s -1 to 4.4 km s _1 , and their spatial distributions are 
shown in Figure 3. The maser spots are distributed over 
an area of 300 mas x 600 mas. Their distribution shows 
three groups, the northern, middle, and southern. The 
red-shifted components with respect to v sys , arc included 
in only the northern and middle groups, and blue-shifted 
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components are concentrated in the southern group. This 
is consistent with the distribution reported by (Motogi 
et al. 2008). And these spots are located in 4." 6 south- 
east of the peak position of sub- millimeter source (SCUBA 
850 fjxa. image was derived from the Canadian Astronomy 
Data Centre), and at the same position within 0."2 as 
the peak of the 8.4 GHz continuum source observed using 
VLA B-configuration (Kurtz et al. 1994). 
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Fig. 2. Scalar averaged cross power spectra of H2O 
maser emission toward IRAS 06061+2151 obtained us- 
ing VERA. The spectra are integrated for 900 s at the 
Mizusawa-Ishigakijima baseline. The dashed line indicates 
the systemic velocity of -1.6 km s _1 (Wouterloot et al. 1989). 



4..1. Measurements of the Annual Parallax 

We used the 1 1 maser spots for a parallax measurement 
(see Sect. 5.1 for details). They were selected as they were 
detected during at least 6 epochs and their size was less 
than 2.5. The least square method was utilized for fitting 
the movement of each maser spot throughout the observ- 
ing epochs to obtain the sinusoidal parallax curve and the 
linear proper motion represented by following equations 
(Green 1985); 



a = Q'o + fi a t + ft Fa 

8 = So + fist + itFs 

where, (ao,#o) is the position at t = 0, ir is the annual 
parallax, (F a , Fg) are sinusoidal functions of direction 
and time describing the parallax ellipse, and (fi a , fis) are 
proper motions in RA and Deck respectively, t is the 
interval in years from the standard epochs. 

We conducted a combined parallax fit, in which the 
positions of 11 spots are fitted simultaneously with one 
common parallax but different proper motions and posi- 
tion offsets for each spot. As the result of this fit, the 
annual parallax of 0.496 ±0.031 mas corresponding to a 
distance of 2.02+o' 12 kp c was obtained (see Table 2, Figure 
4). Besides, the absolute proper motions of 11 maser spots 
are indicated in Figure 3 as black arrows. In the fitting 
process, we obtained the position uncertainties (post-fit 
residuals) of 0.181 mas in RA direction, and 0.183 mas in 
Deck direction, which arc added to make the reduced 
X 2 to be approximately unity. The distance to IRAS 
06061+2151 was determined directly with very high ac- 
curacy of 6.3 %. 

Here we prove that the error bars of 0.181 and 0.183 in 
RA and Deck, respectively are reasonable if we consider 
the uncertainties in astrometry. Possible errors in abso- 
lute position of the target are: a zenith atmospheric delay 
offset at each station, station position error (baseline er- 
ror), and the thermal noise error. A zenith atmospheric 
delay offset at each station is caused by a difference in 
the optical path length (Al) between two sources (target 
source-position reference source) . The range of the eleva- 
tion angle during our observations were between 30° and 
80°. Given the mean zenith angle of observed sources of 
55°, the zenith atmospheric delay residual of 3 cm error, 
and separation of the zenith angles between two sources 
of l.°22, we can estimate Al to be 1.59 mm (Honma et al. 
2007, 2008b; Nakagawa et al. 2008). The positional er- 
ror of 144 ^as due to this effect is derived from Al/B max , 
where -B m ax is the maximum baseline length of VERA. 
The station position error was estimated to be 6 fias for 
a separation angle of l.°22, assuming that the position 
accuracy of 3 mm in each station, based on geodetic ob- 
servations at the S and X bands (Jike et al. 2009). The 
thermal error due to noise was estimated to be 9 fias by 
dividing the half width of the synthesized beam size (9b) 
by the SNR of phase referenced image. We used 0b = 1.8 
mas as the typical synthesized beam size (see Table 1) 
and SNR = 100 as the averaged SNR of phase referenced 
image of maser spots used for astrometric analysis. We 
note that the post-fit residual, which we obtained by the 
least-square fit of annual parallax is consistent with the 
sum of these estimated uncertainties. 

5. Discussion &: conclusions 

5.1. Spot size 

In this section, we discuss the parameter "spot size" 
defined with VEDA (see in Sect. 3). The spot size of 
each maser spot depends not only on the maser structure 
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Right Ascension (J2000.0) 

Fig. 3. 850 (im SCUBA image (gray scale), the yellow cross indicates the center position of 8.4 GHz con- 
tinuum image observed with VLA B-configuration (Kurtz et al. 1994), and the distribution of the H2O 
maser spots. 11 maser spots used for an astromctric analysis arc displayed with open circles and their 
absolute proper motion with arrows. Since they are close one another and they overlaps in this map. 



Table 2. The parallaxes and the absolute proper motions of each maser spot 



Spot ID «lsr Position offset" [mas] Annual parallax b [mas] Absolute proper motion [mas yr 1 ] Detected Epoch 5 N^ pi 







Aa cos 8 


AS 




fi a cos 8 


IH 






1 


-9.480 


-0.325(0.134) 


-0.040(0.132) 


0.518 ±0.084 


0.11±0.12 


-4.04±0.10 


1111*111** 


7 


2 


-9.269 


-0.269(0.197) 


-0.054(0.186) 


0.539 ±0.119 


-0.07 ±0.15 


-3.88±0.12 


1111*1111* 


8 


3 


-9.059 


-0.171(0.190) 


-0.045(0.200) 


0.533 ±0.114 


-0.01 ±0.14 


-3.88±0.13 


1111*1111* 


8 


1 


-8.848 


-0.081(0.173) 


-0.032(0.194) 


0.522 ±0.102 


-0.14 ±0.13 


-3.88±0.13 


111111111* 


9 


5 


-8.637 


-0.033(0.173) 


-0.012(0.191) 


0.495 ±0.103 


-0.19 ±0.13 


-3.89±0.12 


111111111* 


9 


6 


-8.427 


-0.001(0.185) 


-0.005(0.185) 


0.481 ±0.110 


-0.21 ±0.14 


-3.90±0.12 


111111111* 


9 


7 


-8.216 


0.024(0.188) 


0.013(0.169) 


0.473 ±0.111 


-0.20 ±0.14 


-3.92±0.11 


111111111* 


9 


8 


-8.005 


0.049(0.175) 


0.025(0.156) 


0.466 ±0.103 


-0.14±0.13 


-3.92±0.10 


111111111* 


9 


9 


-7.795 


0.050(0.181) 


0.060(0.144) 


0.416 ±0.116 


-0.07 ±0.18 


-4.01±0.12 


111111*1** 


7 


10 


-7.584 


0.114(0.196) 


0.119(0.145) 


0.419 ±0.137 


-0.05 ±0.27 


-3.92±0.18 


111 1 11**** 





11 


-7.373 


18.435(0.129) 


12.334(0.276) 


0.641 ±0.120 


-0.12 ±0.12 


-3.75 ±0.19 




6 


combined fit 




(0.181) 


(0.183) 


0.496 ±0.031 











a: "Position offset" are offsets from the reference position which was set to be a = 06"09 m 06f 970593, <5 = +21°50'41."204897 (J2000.0). We determined this 
coordinate so that the position offset of the strongest maser spot in 2007/296 is (0, 0). Positions are those observed at the epoch detected for the first time, as 
noted in the "Detected Epochs'". Numbers in parentheses represent the uncertainty, 
b: The results were obtained from the RA- Decl. parallax fit. 

c: Each digit shows the detected epochs of the maser spots corresponding to 2007/296, 2007/349, 2008/046, 2008/111, 2008/167, 2008/307, 2008/349, 2009/046, 
2009/131, and 2009/250 respectively. "1" means that the maser spot was detected, and "*" means non-detection, 
d: Total number of detected epochs. 



but on the signal to noise ratio (SNR) of the brightest 
component in the map. Even if the maser spot has no 
structures, low SNR causes an increase of the spot size. 
In order to conduct an analysis of astrometric measure- 
ment with higher accuracy, the maser spots which have no 



structures and whose SNRs are substantially high, should 
be selected. By utilizing the maser spots which falls below 
the threshold for the spot size, it is possible to conduct an 
analysis of astrometric measurement with the maser spots 
which have no structures, or are very bright. 
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Fig. 4. Movements of the 11 maser spots and the results of the least square fit for each maser spot in IRAS 06061+2151. 
Each panel shows as follows; (a) shows the displacement of RA offset against day of year (DOY) from 2007/250, (b) 
shows the displacement of RA offset removed the best fit proper motion from (a), (c) shows the displacement of Decl. 
offset against day of year (DOY) from 2007/250, and (d) shows the displacement of Decl. offset removed the best fit 
proper motion from (c), respectively. To clarify, the data for each spot were added small offset (-0.5 - -32.0 mas) ver- 
tically. The Spot ID of each maser spots is labeled in Table 2 and shown in the most left panel. The error bars 
represent the position uncertainties resulting from systematic errors which are given to make the reduced \ 2 to be ~ 1. 



In this paper, we used the maser spots which were de- 
tected at more than 6 epochs, and whose spot size is less 
than 2.5 for astrometric analysis. On the other hand, 
we also conducted a combined parallax fit for 13 maser 
spots which were also detected at more than 6 epochs, 
and their spot sizes were less than 5.0. As the result of 
the latter analysis, we obtained the annual parallax of 
0.511 ±0.029 mas. Comparing the former result with the 
latter, we found no significant difference. In the analysis 
of this study, there is little difference in the number of 
maser spots which were used for analysis with between 
former and latter conditions, because of the high thresh- 
old set for the number of detected epochs in which the 
maser spots were detected in this analysis. Although var- 
ious values were also used to set the limit for the number 
of detected epochs and the spot size, we could not find 
the result that the structure of the maser spots adversely 
affects astrometric accuracy. The astrometric accuracy in 
present results of this source completely depends on the 
amount of statistics. 

However, if the spot size of maser spot exceeds 3.0, sev- 
eral components appear around the brightest component 
as seen in the most right figure of Figure 1. In the case 
of this map, there would be no effect on the astrometric 
analysis because it is possible to distinguish the brightest 
component from other components (the brightest one ex- 
ists as a compact component). But if some components 
associated with the brightest component are seen in the 
map, such a maser spot would adversely affect the astro- 
metric result and its accuracy. So, we have to confirm 
the shape of maser spots in the phase referenced map to 
determine whether each maser spot has a structure which 



affects the astrometric result, or not. Only if we set a limit 
of 2.5 for the spot size such as in the present study, only 
compact maser spots are used for the astrometric analysis 
as seen in the left and middle maps of Figure 1 . 

5.2. Motion of IRAS 06061+2151 

The systematic motion of maser source is affected by 
its internal structure motions. In order to reduce their 
effects, we first calculated the un-weighted mean of the 
proper motions of the 11 maser spots in Table 2, and 
derived the mean proper motion to be (/Lt a cos<5, fj.$) = 
(-0.10±0.10, -3.91±0.07)masyr- 1 . Here, we calculated 
the standard deviation of the proper motions of 11 maser 
spots as the uncertainty of this mean proper motion. In 
order to obtain the proper motion in the Galactic coor- 
dinate, the tracking center position of IRAS 06061+2151 
and the position added the mean proper motion values to 
the Galactic coordinate. We obtained the proper motion 
with respect to LSR in the Galactic coordinate to be (fxi, 
fib) = (3.37 ±0.01, -1.98 ±0.12) mas yr" 1 by the differ- 
ence of them. Given the source distance of 2.02 kpc, the 
proper motions (/i;, fib) correspond to a velocity of (vi, 
v b ) = (32.20±0.12, -18.92±1.14) km s~\ Also, the sys- 
temic velocity of —1.6 ±0.2 km s _1 derived from 12 CO 
observation was used as the radial velocity of this source. 

In order to convert i>lsr = —1-6 ±0.2 km s~ x to a helio- 
centric frame, i>heik>, we used the following equation with 
the standard solar motion values ([/Q td , Vg td , W@ d ) = 
(10.3, 15.3, 7.7) km s _1 (see the appendix of Rcid et al. 
2009b). 

Whdio = Wlsr - (C/o td cos/ + Vg d sin I) cos b - Wg d sin b 
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*See Sect. 5.2 



As the result of this calculation, we obtained Wheiio to be 
10.8 km s" 1 . 

Next, to discuss about the Galactic 3-dimensional mo- 
tions, we used the solar motion values of (Uq, Vq, Wq) = 
(10.0±0.36, 5.25±0.62, 7.17±0.38) kms^ 1 , derived from 
Hipparcos data by Dehncn & Binney (1998), Rq = 8.5 
kpc as the distance of the Sun to the Galactic center, and 
0o = 220 km s _1 as the rotation speed of the LSR. Also, we 
assume a flat rotation which means the rotation speed at 
the position of the target source &(R) = Qq, where R is the 
distance from the target source to the Galactic center (in 
the case of IRAS 06061+2151, R = 10.5 kpc). As a result 
of calculations using these observed and assumed param- 
eters, we found that IRAS 06061+2151 implies peculiar 
motions of {U p , V p , W p ) = (9.60 + 0.67, -28.50 ±2.26, 
— 11.56+ 1.55) km s _1 . where U p is peculiar motion of a 
star toward the Galactic center, V p is peculiar motion of a 
star in the direction of the Galactic rotation, and W p is pe- 
culiar motion toward the North Galactic Pole. With this 
calculation using the Galactic parameters assumed above, 
IR AS 06061+2151 showed a peculiar motion of larger than 
28 km s" 1 counter to the Galactic rotation. In order to 
confirm whether this peculiar motion was caused by the 
model of the Galactic rotation or not, we re-calculated 
the peculiar motion of IRAS 06061+2151 by using the 
solar motion of (U & , V Q , W Q ) = (10.3, 15.3, 7.7) and 
(ll.lOtg-^, 12.24i°'g, 7.25±° il) kms" 1 (Schonrich et al. 
2010), #0 = 8.0 ±0.5 kpc (Rcid 1993), and 6 = 236±15 
km s _1 (Reid & Brunthaler 2004). We list the result of re- 
calculation of the peculiar velocity components in Table 3. 
Though the Galactic parameters were variously changed, 
the motion of IRAS 06061+2151 counter to the Galactic 
rotation was larger than 17.0 km s _1 . Therefore, we sug- 
gest that IRAS 06061+2151 orbits the Galaxy slower than 
its rotation. 

The peculiar motions of several sources which are lo- 
cated in the Perseus spiral arm are shown in Table 4. 
The parallaxes and the proper motions of these sources 
have already been measured by previous studies. We cal- 
culated (Up, Vp, Wp) of each source using R = 8.5 kpc, 
6 = 220 kms" 1 , {U q ,V q ,Wq) = (10.0 + 0.36, 5.25 + 0.62, 
7.17 + 0.38) km s _1 and the parameters of each source al- 



ready reported (Moellenbrock et al. 2009; Moscadelli et 
al. 2009; Oh et al. 2010; Reid et al. 2009a; Sato et al. 
2008; Xu et al. 2006). We used the systemic velocities 
derived from CO and HCO + observations (Torrefies et al. 
1987; Woutcrloot et al. 1989; Molinari et al. 2002; Lee & 
Jung 2003) as the radial velocity for H2O maser sources, 
because the line width of H2O maser sources is gener- 
ally very broad (from a few tens of km s _1 to 100 km 
s" 1 ). Six sources (IRAS 06061+2151, IRAS 06058+2138, 
AFGL2789, W30H, NGC7538, IRAS00420+5530) orbit 
the Galaxy at least 15 km s _1 slower than the expected 
Galactic rotation, and they move at on average 27 km s _1 
toward the Galactic center and counter to the Galactic 
rotation. 

We note that NGC 281 does not show this peculiar mo- 
tion, this would be because that NGC 281 is slightly offset 
from the Galactic plane and is associated with an expand- 
ing super-bubble (Sato et al. 2008). Thus, the result of 
our observation of IRAS 06061+2151 would be consistent 
with the conclusion that previous studies imply a non- 
circular rotation of the Perseus spiral arm. 
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